Novel thick-target inverse kinematics method for the astrophysical '*C + !?C fusion reaction* 
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The '?C + '°C fusion is one of the most important reactions in modern nuclear astrophysics. The trend and 
magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes. 
However, direct measurement of this reaction is extremely difficult, which makes it necessary to develop indirect 
methods. In this study, the ?*Na + p reaction system was used to study the compound nucleus 2*Mg. We 
employed a thick-target inverse kinematics method combined with the y-charged-particle coincidence technique 
to measure the proton and a exit channels of 74Mg. Technical details of the 7*Na + p thick-target inverse 


kinematics experiment and analysis are presented herein. 
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I. INTRODUCTION 


12C + 12C fusion reactions play a crucial role in vari- 
ous stellar burning scenarios [1-3], particularly during the 
final stages of massive star evolution, Type Ia supernovae 
events [4, 5], and superbursts [6, 7]. Hence, many direct 
measurements, primarily employing charged particles [8- 
10] and gamma ray spectroscopy [11-13], have been con- 
ducted to study nuclear astrophysical reaction rates [14]. 
Because the Coulomb barrier height for the !?C + !?C sys- 
tem, at approximately 7.5 MeV, is significantly higher than 
the Gamow window energy (Ee.m. = 1.5 + 0.3 MeV), the 
cross-section decreases rapidly to below one nanobarn in the 
energy region of interest. This makes it extremely difficult 
to measure the 17C + !?C fusion reaction directly within 
the Gamow energy region. A series of direct measurement 
data [9, 12, 13, 15] suggests the existence of sub-barrier res- 
onances in the '2C + !?C fusion reaction, leading to an in- 
crease in the S-factor as the energy decreases below 2.5 MeV 
in the frame of the center of mass. Therefore, simple low- 
energy extrapolation based on high-energy data cannot accu- 
rately describe the reaction cross-section within the Gamow 
energy region. 

The !?C+!?C fusion reaction rate is mainly contributed by 
12CQ(12C, p)?3Na, 12C(12C, a)?°Ne, and !2C(!2C, n)?3Mg. 
The contribution of the neutron emission channel is marginal 
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because of its negative Q value [16]. In 2018, the Trojan horse 
method (THM) was employed to measure the S' factor of the 
12G + 12C fusion reaction at Eom. < 2.7 MeV [17]. To 
date, this is the only measurement that has entered the Gamow 
window. The astrophysical S* factor derived from the exper- 
iment is much larger than the values of the compilation [18] 
and various phenomenological and microscopic models [19], 
such as wave-packet dynamics (TDWP) [20], and coupled 
channel calculations, such as CC-M3Y+Rep [21, 22]. By 
contrast, Jiang et al. [23] proposed a hindrance model [24, 25] 
that showed the opposite trend in the energy region of inter- 
est. Despite extensive experimental and theoretical efforts, 
the exact behavior of the 12C + !?C fusion reaction remains 
unclear on the existence of resonances or not, particularly in 
the Gamow window. 


Thick-target inverse kinematics (TTIK) is a novel experi- 
mental method that has been widely used in radioactive ion 
beam measurements [26, 27] over the last two decades. In 
the TTIK measurement, the excitation function is obtained in 
a one-shot experiment with a single-beam energy. Despite 
the simplicity of its experimental setup, it has been proven by 
many measurements [28-35] that a satisfactory resolution of 
the excitation function can be obtained owing to the inverse 
kinematics enhancement. 


Proton resonance scattering induced by radioactive sec- 
ondary beams has been investigated [36] with the TTIK 
method at CIAE since 2005. A series of measurements were 
conducted for 12C + p, 13N + p, !"F + p, and ??Na + p using 
stable and radioactive ion beams [37—42]. In this study, we 
extend the conventional TTIK method to complex exit chan- 
nels for the first time using y-charged particle coincidence 
spectroscopy and demonstrate its applicability to simultane- 
ously extracting the excitation functions of different reaction 
channels. 


Il. EXPERIMENT SETUP 


The experiment was performed in the HI-13 tandem ac- 
celerator laboratory [43-47] at the China Institute of Atomic 
Energy (CIAE) in Beijing. Fig. 1 illustrates the experimen- 
tal setup for the 73Na + p thick-target experiments. A beam 
of 110 MeV ?3Na?* ions, with a current of approximately 
0.2 enA, was directed onto a (CH2), target with a thick- 
ness of 5.8 mg/cm?, resulting in an energy loss of approxi- 
mately 66 MeV. Finally, the ?3Na ions were fully stopped in a 
15.7 mg/cm?-thick carbon target. Because a high-energy Na 
beam can easily undergo fusion-evaporation reactions with 
the carbon atoms in the target, a thick carbon target was used 
to measure the background. 


23Na?* beam 


Fig. 1. (Color online) Experimental setup. 


A silicon telescope system was placed at 0° along the 
beamline to measure the protons and a particles emitted from 
the compound nucleus 74Mg. The silicon telescope sys- 
tem consisted of a 70 um double-sided silicon strip detector 
(DSSD), a 1.5 mm multi-guard silicon quadrant (MSQ), and 
a 1 mm MSQ. The DSSD has 16 strips on each side, dividing 
the entire Si surface into 256 pixels. The small detection unit 
of the DSSD provided high-precision angular resolution dur- 
ing the kinematic reconstruction of the charged particles. The 
charged particles of interest originating from the 77Na + p 
reaction primarily deposited energy in the first MSQ, and the 
signals from the second MSQ were used as anti-coincidence 
signals to reduce the background of high-energy protons and 
a particles. The wall of the target chamber was made of or- 
ganic plastic to minimize absorption of the emitted y rays. 
Six 3-inch lanthanum bromide detectors were uniformly ar- 
ranged around the target chamber and placed on sliding rails 
to allow movement along the vertical beam direction to mea- 
sure the characteristic y rays from the residual nuclei 73Na 
and Ne. 


II. y-PARTICLE COINCIDENCE TECHNIQUE 
A. LaBrg detector array calibration 


A three-component y source of ®°Co, !87Cs, and 244 Am 
was used to calibrate the LaBr3 array. The Geant4 pack- 
age [48] was also applied to simulate the efficiency curves 
of the LaBr3 detectors. The simulation of y-ray energy spec- 
trum of °°Co was compared with the measured data, as shown 
in Fig. 2. Using this set of simulation parameters, the ef- 
ficiency of y detectors at different energies was obtained. 
The efficiency curves of the LaBr3 detectors were fitted us- 
ing Eq. (1) [49], considering several points in the simulation, 
as shown in Fig. 3. The systematic error in the simulation 
process was approximately 3 %. 
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Fig. 2. (Color online) Geant4 simulation for ®°°Co y source 
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i=0 
For the residual nuclei 7*Na and ?°Ne, the full-energy peak 


efficiencies at the characteristic energies of 440 and 1634 keV 
were 28.2 % and 12.6 %, respectively. 
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Fig. 3. (Color online) Efficiency curve of the LaBr3 array. The 
black points are the full-energy peak efficiency of the LaBr3 array 
for y rays of different energies obtained from Geant4 simulations. 
The red line represents the fitting result. 


B. --charged particle coincidence 


The y background originates mainly from the fusion- 
evaporation reaction of the high-energy 7?Na beam with a 
carbon target nucleus. For the exit channel pı of the com- 
pound nucleus, its residual nucleus 23Na emits 440 keV char- 
acteristic y rays, that is, corresponding to the first excited 
state of 7°Na. Because of the large number of y-high-energy 
proton accidental coincidence events, it is difficult to quanti- 
tatively analyze the p; exit channel based on purely -y-charged 
particle coincidence. The black and red lines shown in Fig. 
4 represent the y energy spectra associated with the proton 
emissions. After subtracting the C-induced background and 
correcting for a y efficiency of 440 keV, the absolute counts 
of the p; channel were obtained. 


zNa 


10° 


(CH2)n target 
C target 
BNatp 


rrm 


as 
Q 

n 

E440 keV, 


= 
Q 
T 

€—1273 keV, ”Si 


Counts per 6 keV 


= 
O, 
w 


10? 


PT 


ro er gg ge gs ig gl ge gy lg , I 
500 1000 7500 2000 2500 
Ey (keV) 


Fig. 4. (Color online) y ray single energy spectrum measured with 
the LaBr3 detector array in coincidence with protons. The y spec- 
trum obtained with (CH2),-target and carbon-target are presented as 
black and red lines, respectively. The purple line is the net count 
spectrum obtained by subtracting the carbon background spectrum 
after beam normalization. 


IV. TWO-BODY KINEMATICS RECONSTRUCTION 


The energy spectrum of charged particles obtained by a 
silicon detector is a mixture of a series of excitation func- 
tion effects. Because y-charged-particle coincidence intro- 
duces a large number of interfering events originating from 
the carbon-induced background, it is necessary to subtract 
these effects from the energy spectrum. In the energy range 
of the measurement, the exit channels of the compound 24Mg 
consist of po, pı, Qo, and a,. Taking the po exit channel as 
an example, the detailed processing steps for the two-body 
kinematic reconstruction are outlined below. 

Step 1: Extract the energy spectrum of the proton exit chan- 
nels from the total proton energy spectrum obtained by the 
silicon detectors after particle identification and energy cal- 
ibration. The energy spectrum of the pı exit channel can 
be obtained from the charged particle spectrum of the sili- 
con detector via y1-proton coincidence. After subtracting the 
pı energy spectrum corrected by y efficiency from the total 


proton energy spectrum, the energy spectrum of po can be 
obtained. At this stage, both the po and pı exit-channel en- 
ergy spectra contain C-induced backgrounds. Therefore, the 
proton energy spectrum obtained from the carbon target must 
be analyzed in a similar manner. 
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Fig. 5. (Color online) Relationship between FEec.m. and the deposited 
energy in silicon for the pı exit channel. 


Step 2: Perform an event-by-event two-body kinematic re- 
construction based on where the reaction ?3Na(p, py)’ Na 
occurs within the thick target and calculate the relationship 
between FE... and the energy deposited in the silicon detec- 
tors. Through energy-loss calculations and reaction kinemat- 
ics, the energy deposition in the silicon detectors Es; of pro- 
tons is attributed to Ee.m. of the ?3Na(p, pyi)?3Na reaction. 
As shown in Fig. 5, the red and black dots represent the corre- 
spondence between Es; and Ee.m. for the 73Na(p, py)” Na 
reaction on the (CH2)» and carbon target, respectively. 

The correspondence between Es; and Ec.m. can be fitted 
using the least-squares method with a linear function, as ex- 
pressed in Eq. (2). As shown in Fig. 5, the red and black dots 
approximately follow straight lines. 


Ecm. (£) = Esi(x)-k +b, (2) 


Step 3: Subtract the background events from the energy 
spectrum in the frame of the center of mass. Corrections for 
the beam particle and !?C atom numbers are required for the 
carbon-induced background. The (CH2),, and carbon targets 
have different energy stopping powers for 7?Na. The cor- 
rection factor fa for the '2C atom numbers is the ratio of 
the number of C atoms corresponding to 1 keV of energy de- 
posited in the (CH2),, target and the carbon target for the same 
energy of the ?3Na beam. This relationship was calculated 
using LISE++ and fitted to a linear function (Fig. 6. 
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Fig. 6. (Color online) Relationship between the bin target thickness 
and center-of-mass energy. 


Eom. * k(CHa)a as b(CH2)n 
— Fem. n ke T be = 


The number of !?C atoms corresponding to each bin of 
the carbon target and the (CH2), target energy spectra is 
given by Eq. (3) and Eq. (4). The correction factor fn is 
given by Eq. (5). Subsequently, the carbon background can 
be subtracted from the energy spectrum of po. As shown in 
Fig. 7, the black line represents the energy spectrum of the 
(CH), target and the red line represents the energy spec- 
trum of the carbon target. The net energy spectrum of the 
1H(?3Na, p)?3Na channel can be obtained by subtracting the 
normalized carbon background. 
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Fig. 7. (Color online) Yields from 1H(73Na, po)” Na elastic scat- 
tering channel. The background from the pure carbon target is flat 
without any resonance structure. 


Step 4: Calculate the excitation function of 
1H(?3Na, p)?3Na, based on the net yield of po exit 
channel. Because each bin in the pp energy spectrum corre- 
sponds to different H atom numbers of the (CH2),, target, the 
differential cross-section needs to be calculated separately 
for each energy point. The excitation functions of the proton 
and a channels are shown in Fig. 8 and Fig. 9, respectively, 
according to Eq. (6). 


Neount 
Q = ———— 
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where Ny and Ny represent the numbers of H atoms in 
the target and incident 7Na°* beam ions, respectively. Q 
denotes the solid angle of the Si detector. 
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Fig. 8. (Color online) Excitation function for the proton exit 
channels. po and pı represent TH(?3Na, po)” Na and TH(??Na, 
pi)? Naixao, respectively. 
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Fig. 9. (Color online) Excitation function for the @ exit channels. 
ao and a represent TH(?Na, ao)? Ne and 1H(?Na, a1)? Než 634, 
respectively. 


Throughout the data analysis process, each step resulted 
in an error that was ultimately contained in the individual 
data points of the excitation function. These errors mainly 
include ~y-charged-particle coincidence, statistical, carbon 
background deduction, and pı (a1) deduction errors in the 
po (ao) calculations. For y-particle coincidence, the coinci- 
dence efficiency mainly depends on the y detection efficiency, 
silicon efficiency, and accidental coincidence events. The er- 
rors of these factors were independent. The y detection effi- 
ciency error evaluated through the Geant4 simulation process 
was less than 3%, whereas the silicon detection efficiency er- 
ror was negligible. Accidental coincidence events were uni- 
formly distributed in the coincidence time spectrum and were 


evaluated as having an error of less than 2%. Before the car- 
bon background deduction, the energy spectrum must be di- 
vided into individual bins, with statistical errors assigned to 
each bin. For energy points with low counts, the statistical 
error was relatively large, close to 2%. The error after car- 
bon background deduction depends on statistical errors and 
errors introduced during the subtraction process. During car- 
bon background deduction, the ratio of carbon background 
counts to total counts positively correlated with the error in- 
troduced by the subtraction process. As shown in Fig. 7, for 
the po channel, the error after carbon background subtraction 
at 3.8 MeV was 5%. By comparison, for most energy points 
with lower background contributions, the error after back- 
ground subtraction was approximately 3.5%. For the po(ao) 
channel, similar errors caused by data processing were intro- 
duced when subtracting the p;(a,) channel. The y-efficiency 
error must be considered separately when calculating the fi- 
nal error in the p;(q@ 1) channel. Overall, the two subtraction 
processes were the main sources of excitation function errors. 
The total errors for pı and a were less than 8% and 10%, re- 
spectively. po and ao errors were correspondingly larger, that 
is, less than 10% and 17% over the entire energy range. 


The R-matrix theory [50, 51] is a parametric framework 
that describes compound nuclear reactions and theoretically 
describes the resonance phenomena in nuclear reactions. For 
excitation functions that include resonance states, R-matrix 
analysis can decompose the differential cross section into 
three overlapping components: hard-sphere scattering back- 
ground, independent resonances, and interference of multiple 
resonances. This study provided the exit-channel excitation 
functions of the compound nucleus 74Mg populated by the 
?3Na+p entrance channel. The Azure package was used to 
extract a series of exit channel resonance parameters related 
to the '*C+!“C fusion reaction. The fitting results are indi- 
cated by the red lines in Fig. 8 and Fig. 9. The fitting 
results indicate that nearly 50 24Mg resonances are included 
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V. SUMMARY 


The thick-target inverse kinematics method is widely used 
to measure the excitation functions of p(q) elastic and inelas- 
tic scattering induced by radioactive ion beams. In this study, 
the conventional TTIK method was extended to include com- 
plex reaction channels for the first time, which enables the 
simultaneous extraction of the excitation functions of differ- 
ent reaction channels. The high energy resolution and high 
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